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Abstract. We derive and use an extremely simplified formula for the interaction Casimir 
energy for two separate bodies in the weak coupling regime for massless scalar fields. We derive 
(~| ' closed form solutions for a general arrangement of two (5-function plates finite in one direction 

and infinite in another. We examine the situation of two parallel plates finite in both transverse 
directions. 
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' !• Introduction 

Recently, Emig et al. |2] rederived the multiple scattering formalism and used it to calculate 
new results for Casimir energies between disjoint bodies. Similar techniques have been employed 
for many years, perhaps starting with Renne in 1971 [3j, who rederived the Lifshitz formula [4j. 
More recently many papers have used multiple scattering techniques to examine the correction 
. ^ ' to the proximity force approximation for different situations O O [U [9l [TO] . 

^ . Casimir calculations have conventionally been carried out assuming that parallel plates had 

^ I infinite extent. Gies and Klingmiiller have found corrections to the infinite-size approximation 

for the cases of perfectly conducting boundary conditions [H]. We will here use the multiple 
scattering formalism to examine ^-function plate geometries with finite lengths for scalar fields 
subject to very weakly coupled boundaries. 

The weak coupling regime is interesting to study because the formula for the Casimir energy 
simplifies greatly, and in many cases is amenable to closed form solutions. This introduction 
gives a very quick derivation of the weak-coupling form of the Casimir energy for massless scalar 
fields. For a more complete derivation see Kenneth and Klich [12\ I13j . or Milton and Wagner 
[lOj. 

The vacuum expectation value for the action of a scalar field can be given by the standard 
formula 

Vr = -^Trln|-, (1) 

where Q is the Green's function that satisfies the wave equation of the scalar field, including 
any interactions with any objects or potentials. Go is the free Green's function that satisfies the 
wave equation with the same boundary conditions at infinity as the full Green's function but 
without any interaction with background potentials. For a time-independent system we can use 



the condition that W = 



— J dtE to identify the energy as 



^ = A/^TVln^, (2) 

where G is the Fourier-transformed Green's function given by 

Q{x,x') = I (3) 

and the trace if over spatial coordinates. Given two nonoverlapping potentials such that 
V{x) = Vi{x) + V2(x), we can define an interaction energy as the energy of the full system less 
the energy of each potential acting alone, Sjnt = E{Vi + V2) — E{Vi) — E(y2). This simplified 
expression can be written as 

^int = y ^ TV In (1 - G0T1G0T2) , (4) 

where Tj is the scattering matrix for the z*'* potential defined as Tj = — GoVi)~^. 

By formally expanding out the logarithm in (HD, we can think about the formula as describing 
successively more scattering events between the two objects. The weak coupling expansion keeps 
only the first term of the expansion of the logarithm, essentially describing only a single scattering 
between the objects. Additionally, in the weak-coupling regime the scattering matrix can be 
approximated simply by the potential, T ~ y. This results in a very simplified weak-coupling 
single-scattering approximation to the energy {u — > i(^), 

Eint = [ dCTTGoViGoV2. (5) 



2. 2-)-l Spatial Dimensions 

If the potentials are independent of the z direction then we can further simplify the interaction 
energy. By dividing out the infinite length in the z direction we obtain an energy per unit length, 
which we will represent in this paper by the fraktur symbol (E. By Fourier transforming in the 
z direction we get 

where go is given by 

Go{r - r'- ^e''=^(^-^')5o(r-x - r^; k), (7) 
and K is defined by = + fc^. The two dimensional Green's function is explicitly 

go{r_L-r'^;Hc) = ^Ko{K\r±-r'_^\), (8) 
yielding a weak-coupling form for the energy per unit length of 

e: = -^ / ^^d^-/ ^''^ / d^r'Ki{K\r^-r'^\)Vi{rrX)V2{r'^-X). (9) 

In the case that the potentials are independent of the imaginary frequency ( this simplifies even 
further to 

1 f .2MrL)V2{r'^ 



€ = / dV / d^r'^4^^^^^^^. (10) 




Figure 1. Two cylinders of radii a and b, their centers separated by a distance R. 

To demonstrate the simphcity of this formula we will rederive the Casimir energy for two 
cylinders [lOj. Assume that the two cylinders have radii a and b, and their centers are 
separated by a distance R as shown in figure [TJ This situation can be represented by potentials 
Vi = Xi5{r — a) and V2 = \2^{r' — b) where r and r' are radial coordinates in cylindrical polar 
coordinate systems centered on the respective cylinders. Using p^ . these potentials yield an 
energy per unit length 

Wft f'^ 1 

327r3 Jo Jq R"^ + a^ + b'^ -2aR cose + 2bR cos 9' -2abcos{9- 9')' ^ ' 

With a simple change in angular coordinates to u = 9 — 9' and v = this expression can be 
integrated to yield the exact closed form [lOj 

g = -^^-^^"^ , ^ (12) 

StT ^(i?2 _ (a - 6)2) (i?2 _ + 5)2) ^ > 

2.1. General Configuration of Plates 



^ . d 

Figure 2. Two finite plates of length Li and L2 in a general configuration. In a coordinate 
system centered at the edge of the top plate if is the relative angle between the plates, a is the 
perpendicular distance the lower plate is shifted down, and d is the lateral distance the edge of 
the lower plate is shifted to the left. 

Two finite plates of a general configuration as shown in figure [2] can be represented by the 
potentials in cylindrical coordinates, with origin at the left edge of plate 1, 

Vi = Xi6{e - y^) ®^-^^ ~ (13a) 
r 

V2 = X2S{y + a)@{x + d)e{L2 -d-x). (136) 
Here Q{x) is the step function. 




Using these potentials with ([TO]), we get the following expression for the energy per unit length, 



A1A2 
"327r3 
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L2-d 
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[x — r cos ip)"^ + (o + r sin if)"^ 



This integral can be done exactly, yielding a closed form for the general configuration. 
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where Ti2 is the generalized inverse tangent integral defined by 



Ti2 (x, a) 



^ ArcTan y 



dy- 



y + a 



(15) 



(16) 



(17) 



2.2. Torque 




Figure 3. A finite plate of length L, above an infinite plate, The center of the finite plate is 
a perpendicular distance b above the infinite plate, and the finite plate makes an angle ip with 
respect to the infinite plate. 



It would be desirable to examine the rotational stability of parallel plates. This can be done 
simply by looking at the sign of the torque for small angular displacements, while leaving the 
center of mass stationary. Although one can get a general form for the torque from (jl6p . it is 
much simpler to restructure the problem to isolate the torque as shown in figure [3l This is the 
same as the "Casimir pendulum" problem studied by Scardicchio et al., who used the optical 
approximation [15j . 

Given a situation of a tilted plate of length L over an infinite plate (Li — > L, L2 ^ 00, 
d ^ — cxD in figure [2|), we can very easily isolate the torque that the finite plate experiences. The 
energy per unit length is given by 

g=- lnf ^+f""^ V (18) 

327r^ sin 99 \ 5 _ ^ gin / 

The torque per unit length T found by taking the negative derivative of the energy with respect 
to the tilt angle. This gives an expression for the torque on the plates as 

AiA2Cosyp / 1 ^^( b + ^smip \ Lb \ 

327r2 sin(^ \^sin(^ yb-^sm^J 52 _ ^ V / ' 



^ The generalized inverse tangent integral is related to the dilogarithm function, and much information about it 
can be found in [14] , 



From the expression for the torque we can see some clear quahtative features. The torque 
has a zero value at <p = 0, and a quick evaluation shows that the first derivative is positive, of 
value 



dip 



ip=0 



A1A2 L"^ 
192^63 ' 



(20) 



signifying an unstable equilibrium. For values of L such that L > 2b, the torque diverges as 
simp approaches 2b/ L. This is simply the result of the fact that the plates would touch in this 
situation. If L < 26 the torque has another zero at ip = ^, with a first derivative of 



d1 
dp 




(21) 



which is negative for all L < 2b, meaning a stable equilibrium. Therefore a finite flat plate 
suspended above another infinite plate will tend to orient itself perpendicular to the plate if left 
to rotate about its center of mass. Physically, the reason for this is clear: because the magnitude 
of the Casimir force decreases with distance, to minimize the energy for a fixed center of mass, 
the smaller plate wants to rotate so as to place the shortest side closest to the infinite plate. 
(For the situation when the thickness of the plate is finite, see |16j.) 



2.3. Parallel Plates 
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Figure 4. Parallel plates 



Parallel plates are perhaps the most interesting special case. We can compare the exact 
expressions for energy and force to those for infinite parallel plates, getting corrections for finite 
size. In addition the parallel plates case, due to its simplicity, lends itself well to studying both 
normal and lateral forces. 

Consider the same setup as in the general case shown in figure El simply letting ip go to zero, 
as shown in figure figure [H The energy per unit length can be derived directly from (jl6l) by 
using the identity 

lim aTi2(x,a) = / dxArcTanx, (22) 



yielding an integral form for the energy per unit length 
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dxArcTanx. 



(23) 



Although an indefinite integral for the arctangent exists, this form is perhaps more illuminating 
because all the physical quantities are in the limits. The forces, which are given as derivatives 
of the energy, are all given in terms of arctangents. 

Equation (j23p yields closed forms for the normal force between the plates and the lateral force 
experienced by the plates by taking the negative derivative with respect to a or d, respectively. 



The general form of the normal force, defined as = —d<E/da, is 
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(24) 



In the limiting case of the plates getting very close together we expect to recover the result for 
the pressure for infinite parallel plates times the area exposed. By mathematically taking a ^ 0, 
we use the large argument expansion of the inverse tangent, 



ArcTan ix) = — l-o^ + 

2 X 6x-^ 



for X — > 00, 



(25) 



to recover the expected result plus corrections to that result. Because the limiting form of the 
arctangent depends on the sign of the argument, the single general equation can give several 
different answers depending on the size and position of the plates. For the situation shown in 
figure [4] the limiting form is 
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((L2-d) + 0(a3)), 



(26) 



and the first correction is zero. However, if the plates are the same size and aligned the limiting 
form of the force with the first correction is 
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(27) 



If we let one end of both plates extend off into infinity then we can get the edge correction for 
two aligned plates. This correction is 
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(28) 



The general form of the lateral force, similarly defined as = —d(B/dd, is 
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(29) 



From the exact expression for the lateral force, we find there is only one equilibrium position 

' 1-L2 
2 



occurring at d = , where the derivative of the force is negative: 



dd 
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(30) 



signifying a stable equilibrium. The position and qualitative behavior is as expected, the plate 
have an stable equilibrium when they are symmetrically aligned. 



We are also interested in how the lateral force behaves if the plates are very close together. 
To study that we simply take the limit as a ^ 0. Assuming without loss of generality that 
L2 > Li, to lowest order the force is 

+ for d > and > L2 - Li 

g'rf = <( for d > and < d < L2 - Li (31) 



for d < 



This is what we would expect if we approximated the energy simply as the energy per area 
between the two infinite plates times the area exposed between the two plates, and took the 
derivative of this very simple approximation as the force. 

3. Three spatial dimensions 

Until now we have worked in 2+1 dimensions, meaning that the potentials had infinite length in 
one direction. So the finite plates considered in the previous section were more like two infinite 
ribbons. In this section we will work with plates of finite area. 

We start by working out the TGTG formula ([5]) in three dimensions. In three dimensions 
the form of the Green's function is even easier to work with than in two dimensions, 

1 e-lfll^-^'l 

Go(r,r') = -^ -. (32) 

47r \r — r'\ 

This gives a form of the energy 

and if the potentials are independent of Q then the erpression simplifies further to 

If we restrict ourselves to parallel plate 5-function potentials of any general shape the energy 
can written in an even simpler form, 

^1^2 f j2^. f ^2 J r„2 , ^/ ^21-3/2^ ^g^^ 



where Ai and A2 are the areas of the two plates, and a, again, is the separation between the 
plates. If we let one of the areas, let it be ^2, tend to infinity, then the energy for a single finite 
plate above an infinite plate is 

^ = -^^1, (36) 

exactly the energy per area from the Lifshitz formula times the area of the finite plate. For weak 
coupling plates, if one of the plates is infinite, and the other finite, then there is no correction to 
the Lifshitz formula. This is not unexpected, and is a result of the fact that the weak coupling 
approximation is the same as pairwise summation. 
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Figure 5. Two finite parallel rectangular plates. The sizes of the plates are Lix x Liy for 
the top plate, and L2x x for the bottom plate. The lower left corner of the upper plate is 
displaced a distance dx in the x direction and dy in the y. 



3.1. Rectangular Parallel Plates 

For two rectangular parallel plates, as shown in figure EJ the interaction energy is given by the 
integral 
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This expression can be partially integrated and rewritten as 
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The two-dimensional indefinite integral in the equation is given by 



(37) 
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-^xln(l +y2) _ iyin(l + x2) - \/l + x^ + y2. 



(39) 



The final closed-form expression for the energy of the two rectangular parallel plates is somewhat 
messy, consisting of the above indefinite integral evaluated at 16 different combinations of 
variables. 

The normal and lateral forces can again be given by the derivatives of the energy with respect 
to the separation a or to the displacement (this time either dx or dy). 

The lateral force from the plates has a stable equilibrium when the centers of the two plates 
are aligned. However, first derivatives of the force can be different for displacements from the 
equilibrium position in the x and y directions depending on the geometry. 

Perhaps the most interesting property of this system as in section 12.31 to examine what 
happens to the attractive force between the plates as the plates get very close together. For 
very small separations we should get an expression for the force as a power series in a where the 
first term is the pressure given by the Lifshitz formula times the area between the plates. 
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Figure 6. Two coaxial disks, separated by a distance a. The radii for the two disks are i?i and 



Using the large argument expansion for the arctangent (j25p . it is possible to get such an 
expression for the two plate arrangement, although the expressions for the area and the correction 
terms depend on the layout of the plates. For a situation in which the upper plate is completely 
above the lower plate, with none of the edges aligned, the area is given as ^ = Li^Liy and the 



first correction term is ci 
exactly aligned {d^ = dy - 



0. For a situation where both plates are the same size, and they are 



0) then the area is A 
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(41) 



3.2. Parallel Disks 

Instead of asking how two parallel rectangular plates attract, we could just have easily asked 
how two disks attract. The interaction Casimir energy for two coaxial disks separated by a 
distance a is given by 
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(42) 



This expression can be integrated term by term in a series expansion in powers of r and r'. 
The energy can then be expressed as 
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where 



2(m + 1; 
m + 1 
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(44) 



The power series is convergent, so we can simply take the derivative of each term to get the 
force between the plates. By using the asymptotics of the power series in the limit as a — > 
we recover the expected result and get corrections to the infinite plate result. For two different 
sized disks, if Ri < R2 then A = irRf and ci = 0. For two equal sized disks where Ri = R2 = R, 
A = ttR'^ and 

Cl = ^. 45 



4. Conclusions 

The weak-coupling regime greatly simplifies Casimir calculations, and even more so the 
interaction energy between two bodies. The simplicity allows us to obtain closed-form solutions 



to the energy and force between some nontrivial geometries. These closed-form solutions can 
help us in understanding the mechanics of these systems, such as the torque on two finite 
plates studied here. Also these solutions can help us in understanding the limits of certain 
approximations, such as the proximity force approximation or the correction to the attractive 
force between plates of finite size. 

For the three cases of parallel plates studied here, the 2+1 D parallel plates, the rectangular 
parallel plates, and the parallel co-axial disks, we get very similar results in the limit as the 
plates get very close together. If the edges of the plates do not align then the first correction 
term ci = and the area is the area of overlap of the plates. This fact can be reconciled with 
the fact that for a finite plate over an infinite plate the exact result for the attractive force is 
simply the pressure from the Lifshitz formula times the area of the plate. We would expect this 
to be a good approximation if one plate were much larger than the other, which corresponds to 
the first correction being zero. If one plate is even slightly larger, in the limit as the separation 
goes to zero the difference in the size of the plates is still large in comparison to the separation. 

In the case of the edges of the two plates aligning exactly (that is, plates of the same size and 
shape directly above one another) then the area is simply the area of the plates and the first 
correction term takes the form 

1 Perimeter 

= <-"" 

This is a general property of this system, and can be visualized by realizing that in the limit 
as a goes to zero locally along any edge the system will appear to be two semi-infinite plates 
with their edges aligned. Therefore we might expect the correction to be a proportional to the 
perimeter of the plates, with the constant of proportionality given by ()28p . 

In addition, these closed-form results will act as simple test cases for numerical studies. 
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